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the  stochastic  formulation.  In  addition,  this  research  proposes  correction  factors  that  may 
be  needed  when  comparing  dispersion  in  angular  mil  at  multiple  downrange  distances. 
The  significance  of  the  correction  factors  was  demonstrated  through  a  few  application 
examples.  [DOI:  10.1115/1.4003430] 


1  Introduction 

Ballistic  dispersion  is  one  of  the  important  metrics  that  has  been 
used  to  assess  the  performance  of  a  weapon  system.  Over  the  past 
two  decades,  much  research  has  been  conducted  to  characterize 
target  impact  dispersion.  In  general,  a  wide  range  of  ammunitions, 
including  large-,  medium-,  and  small-caliber  projectiles,  along 
with  a  few  different  weapon  systems,  including  air  guns  and  elec¬ 
tromagnetic  guns,  has  been  covered  for  dispersion  study  [1-4]. 
Most  investigations  focused  on  dispersion  analysis  of  experimen¬ 
tal  data  for  a  particular  weapon  system.  Some  published  research 
discussed  ballistic  dispersion  with  the  consideration  of  a  pre¬ 
defined  physical  condition  from  external  environment,  and  others 
addressed  the  phenomenon  with  the  demonstration  of  computer 
modeling  and  simulations  [5,6]. 

Conventionally,  the  absolute  value  of  dispersion  was  initially 
computed  based  on  a  number  of  rounds  at  a  certain  range.  Then, 
the  dispersion  quantity  was  divided  by  the  corresponding  down- 
range  distance  to  obtain  the  dispersion  in  angular  mil.  This  dimen¬ 
sionless  unit  enables  us  to  compare  dispersions  at  two  or  more 
different  downrange  distances.  This  comparison  of  dispersion  at 
multiple  distances  is  valid  under  the  assumption  that  the  projectile 
flies  in  a  vacuum  space.  In  other  words,  no  external  factor  is 
considered  in  the  calculation.  Thus,  one  of  the  objectives  of  this 
study  is  to  formulate  the  random  nature  of  round-to-round  errors 
in  consideration  of  several  varying  environmental  factors  and 
evaluate  if  the  underlying  assumption  for  the  conventional  disper¬ 
sion  calculation  is  appropriate  when  some  external  forces  are 
taken  into  account  during  flight.  This  analytical  approach  links 
dispersion  results  to  projectile  initial  velocity  conditions  or  vice 
versa  from  a  theoretical  point  of  view.  When  one  of  these  two 
measurements  is  available,  the  other  quantity  could  be  derived 
accordingly.  The  proposed  formulas  in  this  study  may  be  utilized 
to  facilitate  performance  assessments  of  weapon  systems. 

This  investigation  begins  with  the  introduction  of  the  most  sim¬ 
plistic  case,  one  where  no  external  force  is  applied  on  a  projectile 
during  flight.  That  is,  the  initial  conditions  of  the  projectile  at  the 
muzzle  remain  the  same  throughout  the  entire  flight  period.  It  has 
been  observed  from  many  experiments  that  the  target  impact 
points  (TIPs)  were  randomly  distributed  even  when  the  same  am¬ 
munition  and  weapon  types  were  fired.  A  reasonable  interpretation 
is  that  the  spread  of  the  TIP  in  the  elevation  direction  is  mainly 
driven  by  the  variation  in  the  initial  velocity  of  the  projectile  in 
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the  elevation  direction  when  it  exits  a  gun  barrel.  In  addition,  the 
spread  of  the  TIP  in  the  azimuth  direction  can  be  assumed  to  be 
primarily  driven  by  the  variation  in  the  initial  velocity  of  the  pro¬ 
jectile  in  the  azimuth  direction.  It  has  been  known  that  the  varia¬ 
tion  in  the  projectile  initial  velocities  may  be  influenced  by  asym¬ 
metric  geometry,  uneven  in-bore  pressure,  imperfect  straightness 
of  gun  tubes,  etc.  The  uncertainty  from  one  round  to  another  can 
be  better  modeled  in  terms  of  probability  distribution  functions. 
When  a  large  number  of  rounds  are  considered,  it  has  been  ob¬ 
served  that  the  initial  velocity  in  each  elevation  and  azimuth  di¬ 
rection  is  approximately  a  Gaussian  distribution.  This  approxima¬ 
tion  of  a  Gaussian  random  variable  to  the  initial  velocity  will  be 
used  throughout  the  study. 

High-speed  projectiles,  such  as  Mach  3  or  higher,  exhibit  vary¬ 
ing  initial  velocities  in  both  elevation  and  azimuth  directions  at 
the  gun  muzzle.  Thus,  the  velocities  are  considered  to  be  random 
variables.  However,  the  initial  downrange  velocity  was  found  to 
possess  a  fairly  low  coefficient  of  variation  [7]  and,  therefore,  can 
be  assumed  to  be  a  constant.  Furthermore,  three  external  factors — 
drag  force,  gravity  drop,  and  crosswind — commonly  considered  to 
influence  dispersion  of  TIP,  are  taken  into  account  in  the  formu¬ 
lation  of  ballistic  dispersion.  It  is  well  known  that  as  the  down- 
range  velocity  decreases  with  the  increase  in  downrange  travel 
distance  due  to  air  drag,  the  flight  time  is  no  longer  linearly  pro¬ 
portional  to  the  travel  distance.  In  addition,  the  effect  of  gravity 
will  accelerate  or  decelerate  the  initial  velocity  in  the  elevation 
direction,  depending  on  whether  the  initial  velocity  is  going 
downward  or  upward.  There  are  also  two  possible  scenarios  when 
considering  crosswind,  i.e.,  the  initial  lateral  velocity  of  a  projec¬ 
tile  being  of  the  same  or  opposite  direction  of  the  crosswind.  With 
the  aid  of  probability  theory  and  statistics,  the  association  of  bal¬ 
listic  dispersion  with  these  physical  effects  is  demonstrated  in  a 
stochastic  approach. 

For  low-speed  projectiles,  experimental  data  have  shown  that 
the  coefficients  of  variation  in  the  initial  downrange  velocity  were 
significantly  higher  than  those  for  high-speed  projectiles.  Thus,  to 
be  more  accurate,  the  initial  downrange  velocity  should  be  treated 
as  a  random  variable.  Some  of  the  known  causes  regarding  the 
phenomena  include  time-varying  frictions  between  obturator  and 
bore  surface,  gun  barrel  manufacturing  tolerance,  granular  shape 
variations  of  propellant  charges,  packaging  deviations  of  each  pro¬ 
pellant  load,  distinctive  flame-spreading  path,  changing  ambient 
temperature,  etc.  Since  the  contribution  of  these  factors  to  the 
variations  in  initial  velocities  is  significant  in  all  three  Cartesian 
directions,  all  the  assumptions  for  high-speed  projectiles  remain 
effective.  As  a  result,  the  initial  downrange  velocity  is  represented 
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by  a  nonzero-mean  Gaussian  random  variable.  The  mean  of  the 
random  variable  is  the  overall  average  of  the  downrange  velocity 
for  the  weapon/ammunition  system. 

Another  major  goal  of  this  research  is  to  propose  correction 
factors  (CFs)  that  may  be  needed  when  comparing  dispersion  in 
angular  mil  at  multiple  downrange  distances.  Application  ex¬ 
amples  of  the  CF  formulation  are  demonstrated.  In  addition,  the 
sensitivity  of  the  CF  to  drag  coefficient,  projectile  mass,  air  den¬ 
sity,  crosswind  speed,  etc.,  is  also  investigated.  A  parametric  study 
is  performed  to  gain  better  understanding  of  the  influencing  pa¬ 
rameters.  It  is  of  interest  to  learn  to  what  degree  these  parameters 
contribute  to  dispersion  and  what  can  be  expected  regarding  the 
properties  of  ballistic  dispersion.  The  CF  is  found  to  be  significant 
in  several  cases  and  therefore  should  be  taken  into  account  when 
dispersion  comparison  is  made. 

Finally,  it  should  be  mentioned  that  with  the  consideration  of 
the  variation  in  the  initial  velocities  in  all  three  Cartesian  compo¬ 
nents,  the  formulation  of  the  dispersion  appears  to  be  quite  com¬ 
plex  but  robust.  This  study  accounts  for  the  effects  of  aerody¬ 
namic  drag,  gravity,  and  crosswind.  It  is  by  no  means  an 
indication  that  ballistic  dispersion  is  influenced  only  by  these  fac¬ 
tors.  For  instance,  the  effects  of  yaw  rate  and  Coriolis  drift  for 
long  range  projectiles  are  ignored.  Caution  shall  be  taken  when 
employing  the  proposed  formulas  for  ballistic  applications  regard¬ 
ing  the  underlying  assumptions. 

2  Conventional  Calculation  of  Ballistic  Dispersion  in 
Angular  Mil 

Traditionally,  when  a  projectile  exits  a  gun  barrel,  the  initial 
downrange  velocity  uZg  is  assumed  to  be  constant,  i.e.,  Var[u.  ] 
=  0,  from  one  round  to  another.  In  addition,  the  drag  force  is 
ignored  such  that  the  downrange  velocity  of  the  projectile  at  any 
time  remains  constant  over  the  entire  flight  period,  i.e.,  vz=vZq.  In 
the  elevation  direction,  the  initial  velocity  vy  is  considered  to  be 
a  zero-mean  Gaussian  variable,  i.e.,  £[i>j,0]=0  and  Var[uV()] 
=  (SD[uVo])2.  In  addition,  the  effect  of  gravity  is  neglected  such 
that  the  elevation  velocity  of  the  projectile  at  any  time  remains  the 
same  before  it  hits  a  target.  Lastly,  the  initial  velocity  vXq  in  the 
azimuth  direction  is  also  considered  to  be  a  zero-mean  Gaussian 
variable,  i.e.,  £[uVo]  =  0  and  Var[u(o]  =  (SD[uA-o])2,  and  the  velocity 
remains  unchanged  when  no  crosswind  is  taken  into  account.  In 
short,  the  net  forces  are  all  zeroes  throughout  the  entire  flight 
period  in  the  downrange,  elevation,  and  azimuth  directions. 

Given  a  downrange  distance  of  dz  ,  the  flight  time  over  the 
distance  can  be  written  as 


Since  dZi  and  vz  are  constants,  t\  is  a  constant  as  well.  Over  the  t1 
period  of  time,  the  travel  distance  of  the  projectile  in  the  azimuth 
direction  can  be  written  as 

dXl  =  vxtx  (2) 

By  taking  the  expected  value  on  both  sides  of  the  equation,  the 
expected  value  of  dXi  can  be  obtained  by 

£[rfri]  =  £[uVo>i  =  0  (3) 

Similarly,  the  variance  of  d: (  can  be  derived  by  taking  the  vari¬ 
ance  on  both  sides  of  the  equation  as  follows: 

Var[dVi]  =  Var[uAo>j  (4) 

It  can  be  easily  understood  that  the  mean  value  of  dXl  is  zero 

because  v,  is  a  zero-mean  Gaussian  variable.  The  variance  of  dr 
xo  xi 

equals  the  variance  of  vXg  multiplied  by  the  square  of  the  flight 

time.  As  a  result,  the  dispersion  in  the  azimuth  direction  at  dZl  can 
be  expressed  as 


d- 

SD[d  ]  =  SD[n  ]f,  =  SD[u  ]—  (5) 

1  0  0  V- 

z0 

The  ratio  SD[dA|]/<fZ|  multiplied  by  a  constant  6400/(2tt-)  is  the 
conversion  of  the  absolute  dispersion  to  angular  mil,  adopted  by 
NATO  countries.  This  process  can  be  repeated  for  another  down- 
range  distance  d7g.  In  this  case,  the  flight  time  can  be  expressed  as 
t2=dzJvz  .  Similarly,  t2  is  a  constant  since  d-2  and  u.  are  both 
constants.  Over  the  t2  period  of  time,  the  travel  distance  of  the 
projectile  in  the  azimuth  direction  can  be  written  as 

dH  =  vxt2  (6) 

The  statistical  properties  of  dx^  can  be  obtained  by  taking  the 
expected  value  and  the  variance  on  both  sides  of  the  equation  as 
follows: 


E[dX2]  =  E[vXo]t2  =  0  (7) 

and 

Var[<42]  =  Var  [vXg\tl  (8) 

As  a  result,  the  dispersion  in  the  azimuth  direction  at  dZg  can  be 
expressed  as 

d, 

SD[d,  ]  =  SD[u  ]f2  =  SD[u  ]— 

vZo 

Combining  Eqs.  (5)  and  (9)  yields  the  following: 

SD  [dXl]  SDK2] 

d '  “  d ' 

Z1  z2 

The  ratios  of  the  x  dispersion  to  the  corresponding  downrange 
distance  are  equal.  This  equation  confirms  the  geometric  relation¬ 
ship,  as  shown  in  Fig.  1.  In  the  elevation  direction,  it  is  also 
straightforward  to  derive  the  relationship  as  follows: 

SD [d  ]  SD[d  ] 

— ^  =  — (u) 

dZl  dZ2 

That  is,  the  ratios  of  the  y  dispersion  to  the  corresponding  down- 
range  distance  are  equal  as  well.  The  aforementioned  geometric 
relationship  also  holds  in  the  elevation  direction.  Equations  (10) 
and  (11)  yield  the  unit  of  dispersion  in  angular  mil  and  have  been 
commonly  used  to  compare  dispersions  at  two  or  more  downrange 
distances. 


(9) 

(10) 


3  Formulation  of  Dispersion  With  a  Constant  vz 

It  has  been  shown  that  for  high-speed  projectiles,  the  coefficient 
of  variation  of  vz  is  as  small  as  0.5%  through  rigorous  studies  of 
225  distinct  barrel  shapes  and  varying  in-bore  pressure  due  to  the 
variations  in  propellant  charges  [7],  Thus,  the  assumption  that 
v-  =  const,  i.e.,  Vaifu-  1  =  0,  should  be  reasonable.  This  section 
considers  the  effects  of  the  external  factors:  drag  force,  gravity 
drop,  and  crosswind.  That  is,  the  flight  velocities,  vz,  vy,  and  vx, 
vary  over  time.  In  other  words,  the  significance  of  the  three  time- 
varying  in-flight  velocity  components  is  taken  into  account  in  the 
dispersion  formulation.  The  initial  velocities  in  the  azimuth  and 
elevation  directions  are  represented  by  zero-mean  Gaussian  ran¬ 
dom  variables,  i.e.,  £[uV()]  =  0  and  Var[u},o]  =  (SD[uVo])2,  and 
£[ur()]  =  0  and  Var[uto]  =  (SD[i;Ao])2,  respectively,  throughout  this 
study. 

3.1  Downrange  Direction.  The  downrange  velocity  of  flying 
projectiles  decreases  over  time  because  of  air  drag  force.  As  a 
result,  the  ratio  of  d\  to  t\  is  no  longer  equal  to  that  of  d2  to  t2.  It 
has  been  shown  that  the  aerodynamic  drag  force  can  be  written  as 
[8] 
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Scatter  of  projectile  impact  points 
at  4,  range  in  azimuth  direction 


F=-\pACdV2  (12) 

where  p=air  density,  A=reference  area,  Crf=drag  coefficient,  and 
V2=v2  +  v2,  +  u2= launch  velocity  of  the  projectile.  Since  vz^>vy 
and  vz^vx,  the  equation  in  the  downrange  direction  is  expressed 
as 


E[Vy\  =  ■ 


£K]+4:(1“e2CX) 


and 


gvz 

2CK 


(l-e2c>z) 


(19) 


d}h 

dt 


=  -c>: 


(13) 


in  which  Cd=pACdl2m.  The  parameter  Cd  is  a  constant  for  the 
same  ammunition  provided  that  the  shape  and  the  mass  of  the 
projectile  are  all  alike  from  one  round  to  another.  Subsequently, 
the  velocity  can  be  derived  as 


1  +  Ov 


(14) 


As  a  result,  the  downrange  travel  distance  can  be  expressed  as 


SD[u,,]  =  — SD[uv  ]  (20) 

V-  '  “ 

zo 

The  expected  velocity  increase  due  to  gravity  is  shown  in  Eq. 
(19).  In  addition,  the  velocity  path  is  a  straight  line  along  the 
range.  That  is,  the  ratio  of  the  elevation  velocity  to  the  downrange 
velocity  is  independent  of  flight  time.  Subsequently,  the  travel 
distance  in  the  elevation  direction  can  be  expressed  as 


d=dA  -^  +  - 


2CX 


g(l-e2c^)\ 

4 Cfvi  ) 


(21) 


ln(l  +  C*dvzt) 

dz= - ~  (15) 

The  flight  time  t  can  be  obtained  by  rearranging  Eq.  (15)  as 

t=~^(ec>z- 1)  (16) 

Cdvz  o 

The  flight  time  t  appears  to  be  nonlinearly  proportional  to  the 
distance.  However,  t  is  not  a  random  variable.  For  projectiles  with 
a  constant  vz  ,  the  flight  time  t  is  determined  given  a  downrange 
distance,  dz. 

3.2  Elevation  Direction.  When  a  projectile  exits  a  gun  tube, 
it  generally  possesses  an  initial  velocity  in  the  elevation  direction. 
This  velocity  is  assumed  to  be  a  zero-mean  Gaussian  random 
variable.  The  variation  in  the  initial  velocity,  vyQ,  from  one  round 
to  another  is  considered  to  cause  the  spread  of  the  TIP  at  any 
range.  The  scenario  that  the  initial  velocity  goes  downward,  i.e., 
toward  the  earth,  is  discussed  herein.  Based  on  Eq.  (12)  and  the 
inclusion  of  the  gravity  effect,  the  equilibrium  equation  in  the 
elevation  direction  can  be  expressed  as 

C-^  =  -Cdvzvy-g  (17) 


The  flight  velocity  of  a  projectile  in  this  case  can  be  solved  as 


2CX 


(l_e2cX)| 


(18) 


By  taking  the  expected  value  and  the  variance  on  both  sides  of  the 
equation,  it  yields 


Again,  by  taking  the  expected  value  and  the  variance  on  both 
sides  of  the  equation,  the  statistical  properties  of  dy  can  be  derived 
as 


E[dy\  =  dzgl(2cyz)  +  g(  1  -  e2c>z)l(4Cd2 v2o)  (22) 

Var[dv]  =  Var[uVo]|^  j  (23) 


SD[rfv]  =  SD[uV()]—  (24) 

Vzo 

The  expected  projectile  drop  distance  from  the  aim  point  due  to 
gravity  is  shown  in  Eq.  (22).  Equation  (24)  indicates  that  the 
dispersion  in  the  elevation  direction  is  the  standard  deviation  of 
the  initial  elevation  velocity  multiplied  by  the  ratio  of  the  down- 
range  travel  distance  to  the  initial  downrange  velocity.  In  other 
words,  the  dispersion  is  still  linearly  proportional  to  the  down- 
range  distance.  Readers  can  find  that  the  result  here  is  the  same 
when  compared  with  that  in  conventional  angular  dispersion  cal¬ 
culation,  as  shown  in  Eq.  (9).  Thus,  air  drag  and  gravity  are  con¬ 
cluded  to  have  no  influence  on  dispersion  as  vZg  is  unchanged 
from  one  round  to  another.  Same  conclusion  can  be  drawn  when 
the  initial  velocity  is  in  the  opposite  direction  of  the  gravity.  One 
can  easily  derive  the  equations  for  the  scenario.  The  relationship 
between  the  downrange  distance  and  the  expected  drop  distance 
due  to  gravity,  along  with  the  distribution  of  the  TIP,  is  illustrated 
in  Fig.  2.  The  mean  point  of  impact  (MPI)  represents  the  expected 
impact  location  for  the  group  of  rounds. 


Journal  of  Applied  Mechanics 


MAY  2011,  Vol.  78  /  031015-3 


Downloaded  15  Oct  2012  to  128.63.78.81 .  Redistribution  subject  to  ASME  license  or  copyright;  see  http://www.asme.org/terms/Terms_Use.cfm 


d 


d 


Fig.  2  Illustration  of  projectile  drop  distance  due  to  gravity  as 
the  initial  velocity  going  downward 


3.3  Azimuth  Direction.  When  a  projectile  is  fired  in  a  cross- 
wind,  its  velocity  relative  to  the  moving  air  mass  has  a  small 
component  in  the  cross-range  direction,  which  causes  the  projec¬ 
tile  to  deflect  in  the  direction  of  the  wind.  Based  on  Eq.  (12)  and 
the  inclusion  of  crosswind  effect,  the  equilibrium  equation  in  the 
azimuth  direction  can  be  expressed  as 

^7 l  =  -CjVz(vx-wx)  (25) 

at 


Similarly,  the  initial  velocity  of  a  projectile  in  the  azimuth  direc¬ 
tion  could  be  in  the  same  or  opposite  direction  of  the  crosswind. 
This  section  discusses  the  scenario  when  wx  and  vx  in  the  same 
direction  only.  Readers  can  easily  derive  the  equations  as  they  are 
in  the  opposite  direction.  The  flight  velocity  in  the  azimuth  direc¬ 
tion  can  be  expressed  as 


:  — vr,  +  wr\  1  -  — 


(26) 


It  can  be  seen  that  this  equation  shows  no  cross-range  velocity  due 
to  crosswind  initially  when  vX(]=0  and  vz=v-Q.  As  vz  decreases 
downrange,  the  wind  effect  increases  even  when  the  wind  speed 
stays  the  same.  In  most  cases,  the  wind  speed  on  flying  projectiles 
varies  from  one  round  to  another.  However,  when  a  large  number 
of  rounds  are  being  studied,  the  assumption  of  a  zero-mean 
Gaussian  distribution  of  the  crosswind  should  be  reasonable.  That 
is,  the  crosswind  could  be  stronger  on  some  rounds  and  weaker  on 
some  other  rounds.  Nevertheless,  due  to  very  short  flight  periods, 
the  wind  speed  for  each  individual  round  may  be  considered  as  a 
constant  for  that  particular  round.  By  taking  the  expected  value 
and  the  variance  on  both  sides  of  Eq.  (26),  it  yields 


E[v  J  =  E[v  ,1  -*■  +  E[wx]  1  ss*  —  =  0  (27) 


and 


VarfuJ  =  Varft;*  ]  —  I  +Var[wJ  1-—  +2  1 


v.  \  v. 

■  —  — Cov[u_  ,wA] 


(28) 


When  no  or  low  correlation  between  the  initial  x  velocity  and  the 
crosswind  speed  is  expected,  the  covariance  term  may  be  dropped 
from  the  equation.  The  travel  distance  in  the  azimuth  direction  can 
then  be  expressed  as 


d,  ec<idz  -  C*jdz  -  1 

dx  =  ~^vx  + - — — - wx 

Cdvzn 


(29) 


Similarly,  by  taking  the  expected  value  and  the  variance  on  both 
sides  of  the  equation,  it  yields 


Fig.  3  Illustration  of  the  increase  in  the  dispersion  of  target 
impact  points  due  to  crosswind 


d-  r  n  ec^-Cld,- 1 

E[dx]  =  —E[vxA  + - — — - £R]  =  0  (30) 


d  V  I  ec‘<dz  -  C*,d  -  1 

Var[r/J  =  Var[uA  ](  —  I  +  Var[vvJ 


Cdvz  o 


ec^-C*dz-l 


Covk,wJ 


(31) 


Equation  (30)  implies  that  crosswind  has  no  influence  on  the  ex¬ 
pected  distance  deviation.  However,  it  increases  the  variance  by 
two  additional  components,  as  shown  in  Eq.  (31).  Since  the  cova¬ 
riance  term  is  usually  small  and  therefore  neglected,  the  disper¬ 
sion  in  the  azimuth  direction  can  be  written  as 


SDK]  - 


+  Var[wx] 


-  Cjdz  -  1 


2 

(32) 


Equivalent  to  Eq.  (29),  the  travel  distance  in  the  azimuth  direction 
at  any  time  t  can  also  be  expressed  as  dx=vXfdz/vz  +wx(t 
- dzlvz ).  When  an  indoor  experiment  takes  place,  i.e.,  no  cross- 
wind  is  in  effect  for  a  certain  period,  a  time  shift  may  be  simply 
appended  in  the  second  component  of  the  equation. 

Regardless  of  the  relative  direction  between  the  initial  lateral 
velocity  and  the  crosswind,  the  drift  distance  caused  by  the  cross- 
wind  is  zero.  The  mean  of  dx  is  not  affected  by  the  crosswind 
given  the  assumption  that  wind  speed  is  normally  distributed  and 
has  a  zero  mean  across  a  large  number  of  rounds.  In  addition,  it  is 
certain  that  the  crosswind  effect  will  always  yield  greater  disper¬ 
sion  value,  as  shown  in  Eq.  (31).  Figure  3  illustrates  the  un¬ 
changed  MPI  at  any  downrange  distance.  The  new  distributions  of 
the  impact  points  due  to  crosswind  are  displayed  in  a  solid  curve 
along  with  the  distributions  in  a  dashed  curve  when  crosswind  is 
neglected. 

Based  on  Eq.  (32),  the  dispersions  in  angular  mil  at  any  two 
downrange  distances,  dz  and  dz^  can  be  shown  as  Eq.  (33).  In  a 
nutshell,  given  a  certain  range,  the  dispersion  in  the  azimuth  di¬ 
rection  can  be  determined  by  the  variance  of  the  initial  azimuth 
velocity  of  the  projectile,  the  variance  of  the  crosswind  velocity, 
and  the  initial  downrange  velocity, 
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SD  [dXl] 


,  dz  \z  /ecA-CX,  -  1 

Var[yJ  —  +Var[wJ 


CdPza 


SD[42] 


d7.V  /ecA-CX-l 


Var[uJ(  -1  )  +  Var[wJ 


CX, 


(33) 


4  Formulation  of  Dispersion  With  a  Varying  r 

It  has  been  observed  that  projectile  launch  velocity  in  the  down- 
range  direction  actually  deviates  from  one  round  to  another.  This 
is  particularly  true  for  lower  speed  projectiles  that  exhibit  a  higher 
coefficient  of  variation  from  experimental  data.  As  a  result,  the 

case  that  v-  is  considered  as  a  random  variable  should  be  ad- 

*■0 

dressed.  Due  to  the  stochastic  nature  of  vZg ,  at  a  certain  range,  the 
varying  uZg  would  lead  to  uncertain  flight  time,  which  in  turn 
affects  the  outcome  of  the  dispersion  in  the  elevation  and  azimuth 
directions.  This  section  intends  to  formulate  the  dispersion  in  a 
more  comprehensive  situation. 

4.1  Downrange  Direction.  Owing  to  the  variation  of  vz  ,  the 
flight  time  t  at  a  given  distance  will  be  a  random  variable.  As  a 
result,  the  statistics  of  t  can  be  obtained  by  taking  the  expected 
value  and  the  variance  on  both  sides  of  the  equation.  It  yields 

ec*A-  1 

E[t]  =  ~ — £[<]  (34) 

'-'d 

and 

eCddz  -  1  , 

SD[f]  =  — —SD[u-zlo]  (35) 

w 

The  first  equation  implies  that  the  expected  value  of  the  flight  time 
is  proportional  to  the  mean  of  the  inverse  of  the  initial  downrange 
velocity.  This  makes  sense  since  faster  projectiles  require  shorter 
flight  time  to  reach  a  target  at  the  same  distance.  In  addition,  the 
variance  of  the  flight  time  is  proportional  to  the  variance  of  the 
inverse  of  the  initial  downrange  velocity. 

4.2  Elevation  Direction.  When  v .  is  considered  a  constant, 

zo 

gravity  and  air  drag  are  found  to  have  no  effect  on  dispersion.  As 

v-  varies  from  one  round  to  another,  the  variation  results  in  a 

zo 

random  effect  of  gravity  drop  on  a  projectile  at  a  given  range.  In 
other  words,  it  is  a  completely  stochastic  process  over  the  projec¬ 
tile  trajectory.  Again,  this  section  discusses  the  scenario  of  initial 
elevation  velocity  going  downward.  The  case  as  the  initial  veloc¬ 
ity  going  upward  is  left  to  readers. 

The  statistical  properties  of  the  y  velocity  can  be  obtained  from 
Eq.  (18)  for  this  case.  Since  vy  and  v-g  are  both  random  variables, 
it  yields 

E\vy\  =  +  “  e2CA)£K02])  (36) 

Var[uT]  =  v2(  Var^X1]  +  -^(1  -  e^Var^-2]  +  JL(l 

\  d  C d 

-  e2cY)Co v[uTou-\u:2]  j  (37) 

Due  to  the  variation  in  the  initial  downrange  velocity,  the  nonzero 
expected  value  of  the  elevation  velocity  has  an  additional  compo¬ 
nent,  which  is  the  correlation  between  the  initial  elevation  velocity 
and  the  inverse  of  the  initial  downrange  velocity.  The  value  may 


be  significant  from  the  physical  point  of  view.  Also,  based  on  Eq. 
(21),  the  statistics  of  the  travel  distance  in  the  elevation  direction 
can  be  obtained  by 


Rdy  1  =  d-E[vyovZo  ]  +  — ;  \dz  + 


1  -  e2CX 

2c: 


E[vf]  (38) 


where 


and 


1*1  +  — 
4Q 

g  (  j  1  -  e2c^’~ 

2  C* 


i  _  y-ca 


+  dzc*  d~  + 


v%1] 


Uv  V- 1,v2] 
ro  zo  zo J 


‘'V'V 


^z„4] 


Vy  V  / 

>o  zo- 


l*+  2 c:  ) Var[y 

CHtVzo’<] 

(39) 

!]-ok<])2 

(40) 

■(^K2])2 

(41) 

(42) 

It  can  be  observed  that  due  to  the  variation  in  v-  .  the  distribution 

zo’ 

of  the  impact  points  no  longer  exhibits  a  zero  mean,  as  indicated 
in  Eq.  (38).  The  expected  deviation  from  the  aim  point  can  be 
determined  accordingly.  In  addition,  the  quantities  VarfD^X] 
and  Var[i>“2]  are  essential  for  the  variance  calculations  of  projec¬ 
tile  velocity  and  travel  distance.  As  mentioned,  the  covariance 
component  is  usually  small  and,  thus,  neglected.  As  a  result,  the 
formulation  of  dispersion  can  be  approximated  by 


SDK,]  “ 


d2y  wivy ..IE1]  +  7^2  (  d. 


4  c: 


1  —  e2c<^z 

2  c: 


2 

VarK2] 


(43) 


Note  that  a  sign  difference  may  be  obtained  in  the  first  component 
of  Eq.  (38)  as  the  initial  velocity  going  upward,  indicating  dis¬ 
similar  expected  gravity  drop  between  the  cases.  Although  the 
sign  difference  also  takes  place  in  the  covariance  term  of  Eq.  (39), 
the  approximation  of  dispersion  shown  in  Eq.  (43)  is  identical  as 
the  covariance  component  is  neglected. 

4.3  Azimuth  Direction.  When  the  initial  azimuth  velocity  of 
a  projectile  is  in  the  same  direction  of  crosswind,  the  statistics  of 
the  flight  velocity  can  be  derived  from  Eq.  (26).  It  yields 


E[v. J  =  vz{E[vXqv^\  +  E[  wxvZg])  (44) 


Var[uJ  =  ^(VarK^"1]  +  Varfny;"1])  +  Var[wJ 

-  2v2zCo\\vxjv~^,wxv~^\  +  2uz(CovKou~1,wJ 

-  Covfw^.wJ)  -  v2zCov\vXov-^,wx,wxv~zl\  (45) 

It  is  shown  that  due  to  changing  vZg  from  one  round  to  another,  the 
expected  velocity  in  the  azimuth  direction  is  no  longer  zero.  Be¬ 
cause  of  three  variables  being  considered,  the  variance  of  the  azi¬ 
muth  velocity  is  fairly  complex.  However,  for  simplicity,  one  may 
ignore  the  covariance  terms  since  the  values  are  usually  small 
compared  with  the  variance  quantity.  Subsequently,  the  statistics 
of  the  travel  distance  in  the  azimuth  direction  can  be  obtained 
from  Eq.  (29).  It  yields 


ecAz  -  C*A  -  1 

£[rfj  =  d7E[vxv~']  + - (46) 

6./ 
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Ww[dx\  =  d1?Wc\v\vxvz')\  + 


ecX-cX-l 


Var[ny;  '] 


+  2  d. 


ec*A  -  C*jdz  -  1 

C 


Cov[u  ,  Wj.11^1]  (47) 


where 


Co^V’v  "'A)1]  =  Elvx0w*v$  ~  (48) 

As  expected,  the  nonzero  mean  of  the  velocity  leads  to  the  non¬ 
zero  mean  of  the  drift  distance,  as  shown  in  Eq.  (46).  When  the 
covariance  term  of  Eq.  (47)  is  neglected,  the  dispersion  can  be 
expressed  as 


SD[rfJ=  \/rf^ar[tV;1]  + 


ecA-C^z-  1 


C, 


Vaifw^1] 


(49) 


As  the  initial  azimuth  velocity  and  the  crosswind  are  in  the  oppo¬ 
site  directions,  the  results  were  found  to  be  similar  to  those  in  the 
elevation  direction  as  discussed  previously. 


5  Correction  Factors  for  Dispersion  Comparison  at 
Multiple  Ranges 

5.1  Projectile  With  a  Constant  pz  .  It  has  been  shown  that 
with  a  constant  u.  ,  dispersion  in  the  elevation  direction  is  linearly 
proportional  to  the  downrange  distance  even  when  air  drag  and 
gravity  are  considered.  In  this  case,  the  geometric  relationship 
shown  in  Fig.  1  still  holds,  indicating  that  the  conventional  ap¬ 
proach  to  obtaining  angular  dispersion  and  making  comparisons  at 
multiple  ranges  is  appropriate.  However,  when  crosswind  is  con¬ 
sidered,  the  dispersion  at  a  longer  range  tends  to  be  overly  esti¬ 
mated  in  the  azimuth  direction.  Thus,  a  CF  when  comparing  dis¬ 
persion  at  two  downrange  distances,  dz  and  dzy  is  proposed  as 


SDK]  SDK], 


-(CF) 


(50) 


Based  on  Eq.  (33),  the  relationship  can  be  written  as 


SD  fe,]  SD  [rf,2], 
d '  d . 


Var[irj(  ^  )  +Var[wJ 


-  cX,  - 1 


CX 


d7 

Var[uJ(^)  +Var[wJ 


ecA-CX2~  1 


2d„ 


(51) 


where  the  components  in  the  curly  bracket  represent  the  CF.  It  can 
be  observed  that  the  CF  depends  on  Var[ut()],  Var[wJ,  and  Cd. 
Given  any  two  downrange  distances,  the  CF  can  be  calculated 
accordingly. 

This  section  adopts  a  5.56  mm  NATO  ammunition  for  paramet¬ 
ric  study  since  it  has  been  under  firing  test  at  the  Army  Research 
Laboratory  over  the  past  several  years.  The  5.56  mm  ammunition 
is  a  current  U.S.  Army  standard  issue  and  is  manufactured  in 
significantly  large  quantities.  Some  of  the  nominal  properties  re¬ 
quired  for  dispersion  calculation  are  given  in  Table  1.  In  light  of 
the  information,  the  parameter  Cd  is  calculated  to  be  0.001151. 
The  absolute  dispersions  in  the  elevation  and  azimuth  directions  at 
200  m  and  500  m  ranges  are  available  from  previous  analysis  and 
are  provided  in  Table  2  [9].  It  should  be  mentioned  that  the  5.56 
mm  firing  test  was  conducted  in  an  indoor  facility  having  a  length 
of  300  m  where  no  crosswind  was  considered. 

From  a  previous  analysis,  the  traditional  approach  yielded  dis¬ 
persions  of  0.2526  and  0.3630  (in  angular  mil)  at  200  m  and  500 
m,  respectively,  in  the  elevation  direction.  If  the  initial  downrange 
velocity  is  assumed  to  be  unchanged  from  one  round  to  another, 
no  correction  is  needed.  Additionally,  the  traditional  approach 
yielded  dispersions  of  0.2949  and  0.3559  (in  angular  mil)  at  200 
m  and  500  m,  respectively,  in  the  azimuth  direction.  The  variances 


Table  1  Nominal  properties  of  M855  ammunition 


Launch  downrange  velocity  (m/s) 

927 

Gravitational  acceleration  (m/  s2) 

9.81 

Bullet  reference  area  (m2) 

2.55X10- 

Mass  (kg) 

0.004 

Air  density  (kg/m3) 

1.204 

Drag  coefficient 

0.3 

of  the  initial  azimuth  velocity  and  the  crosswind  speed  could  be 
back-calculated  by  Eq.  (32),  and  they  were  determined  to  be  0.053 
and  0.23,  respectively.  As  a  result,  a  correction  factor  of  0.732  is 
needed  such  that  the  comparable  dispersion  should  be  0.2949  and 
0.2605,  respectively.  When  a  low  coefficient  of  variation  of  1%  is 
used  for  vz  ,  through  Monte  Carlo  simulation,  a  correction  factor 
of  0.969  is  needed  in  the  elevation  direction  such  that  the  compa¬ 
rable  dispersion  should  be  0.2526  and  0.3518. 

In  addition,  it  is  also  of  interest  to  assess  the  degrees  of  contri¬ 
bution  to  correction  factors  by  parameters  such  as  Cd,  VarfiqJ, 
and  Var[wJ.  The  C*d  mainly  depends  on  drag  coefficient,  bullet 
mass,  and  air  density.  Thus,  possible  values  of  each  factor  are 
adopted  to  determine  C*d,  as  shown  in  Table  3.  The  calculated  Cd 
values  range  from  0.001151  to  0.001919.  It  can  be  observed  that  a 
higher  drag  coefficient  and  air  density  yield  higher  Cd.  In  addition, 
a  lighter  projectile  mass  yields  a  greater  C*d  value.  Overall,  the 
changes  in  drag  coefficient  are  predominant. 

A  total  of  seven  cases  with  different  parameter  values  are 
adopted  and  shown  in  Table  4.  The  CF  at  five  different  ranges,  i.e., 
200  m,  300  m,  400  m,  500  m,  and  600  m,  are  determined  against 
the  baseline  dispersion  at  100  m.  Table  5  gives  the  CFs  for  the  five 


Table  2  Standard  deviation  of  M855  impact  points  about  MPI 
at  200  m  and  500  m  ranges 


Transverse  directions  at  target  location  Standard  deviation 


Elevation  distance  (m)  at  200  m  0.0579 

Elevation  distance  (m)  at  500  m  0.1747 

Azimuth  distance  (m)  at  200  m  0.0496 

Azimuth  distance  (m)  at  500  m  0.1782 
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Table  3  List  of  possible  values 
projectiles 

of  Cd 

for  high-speed 

Case 

cd 

m 

P 

c-d 

i 

0.4 

0.004 

1.204 

0.001535 

2 

0.3 

0.004 

1.204 

0.001151 

3 

0.5 

0.004 

1.204 

0.001919 

4 

0.4 

0.0038 

1.204 

0.001616 

5 

0.4 

0.0042 

1.204 

0.001462 

6 

0.4 

0.004 

1.146 

0.001461 

7 

0.4 

0.004 

1.342 

0.001711 

Table  5  Correction  factors  against  dispersion  at  100  m  in  azi¬ 
muth  direction 


Case 

200  m 

300  m 

400  m 

500  m 

600  m 

i 

0.989 

0.969 

0.938 

0.896 

0.844 

2 

0.994 

0.983 

0.967 

0.945 

0.917 

3 

0.982 

0.948 

0.898 

0.832 

0.755 

4 

0.909 

0.785 

0.659 

0.549 

0.456 

5 

0.999 

0.997 

0.993 

0.988 

0.980 

6 

0.999 

0.997 

0.993 

0.988 

0.980 

7 

0.651 

0.446 

0.327 

0.250 

0.198 

Table  4  List  of  parameter  values  for  the  case  of  constant  vZo 


Case 

c. ; 

Var[uJ 

Var[wJ 

i 

0.001535 

1 

1 

2 

0.001151 

1 

1 

3 

0.001919 

1 

1 

4 

0.001535 

0.1 

1 

5 

0.001535 

10 

1 

6 

0.001535 

1 

0.1 

7 

0.001535 

0.1 

10 

ranges  in  the  azimuth  direction.  Note  that  when  one  is  comparing 
dispersions  at  200  m  and  500  m,  the  ratio  of  the  two  CFs  from  the 
table  is  the  calibration  value  that  should  be  used. 

A  lower  CF  value  implies  that  the  nonlinear  effects  of  drag, 
gravity,  and  crosswind  are  more  significant.  A  CF  near  1.0  sup¬ 
ports  the  usage  of  the  traditional  approach.  In  Table  5,  it  can  be 
seen  that  a  higher  drag  coefficient,  a  lighter  bullet,  or  dense  air 
each  decreases  CF.  The  variance  of  crosswind  speed  appears  to 
influence  CF  significantly,  as  shown  in  case  7  where  CF  deterio¬ 


rates  dramatically.  Overall,  the  correction  factors  decrease  with 
the  ratio  of  the  variances  of  crosswind  to  the  initial  azimuth  ve¬ 
locity.  When  the  two  variances  proportionally  increase  or  de¬ 
crease,  the  correction  factors  remain  the  same,  as  seen  from  com¬ 
paring  cases  5  and  6.  Not  surprisingly,  when  the  variance  of 
crosswind  stays  at  the  same  level,  the  increase  in  the  variance  of 
initial  azimuth  velocity  improves  the  correction  factors.  In  fact,  a 
CF  of  1.0  is  anticipated  when  no  crosswind  is  in  effect. 

5.2  Projectile  With  a  Varying  vz  .  It  has  been  shown  that 
with  a  varying  vZq,  the  angular  dispersion  at  a  longer  range  tends 
to  be  overestimated  in  both  elevation  and  azimuth  directions. 
Therefore,  CF  will  be  formulated  in  both  cases.  In  the  elevation 
direction,  one  can  write  a  comparison  equation  such  that 


SD[dL]  SD [d  ] 

— ^  =  — HCF) 


(52) 


Based  on  Eq.  (43),  given  two  different  downrange  distances,  the 
relationship  can  be  written  as 


SD  [dyi]  SDK2] 


d7 


2  j  _  g2 CjdZi  \  1 

v"[o?k 


l-e2CA 


2  C*, 


Var[u-o2] 


(53) 


The  expression  in  the  curly  bracket  represents  the  CF.  It  can  be  seen  that  the  quantities  Varfu^u^1]  and  Var[u.“]  need  to  be  calculated 
in  order  to  determine  the  CF.  Similarly,  in  the  azimuth  direction,  based  on  Eq.  (49),  the  relationship  at  two  different  downrange 
distances  can  be  obtained  by 

r  x 


SD[rf  ]  SD[d  ] 
- —  = - — 

d7  d, 

Z1  z2 


dzVar[vXov;o]  + 


;CA  -  C*,d-  -  1 


VarKu-1]^ 


d7  Var[vxv;(i  ]  + 


.  eLA-CX  -  1 

r2  \T _ r  -l"1  I  d  Zr> 


2  dZi 

VarjwyC1] 


(54) 


The  expression  in  the  curly  brackets  represents  the  CF.  It  can  be 
seen  that  the  quantities  Var[uV()u7j]  and  VarfWjUj/]  need  to  be 
calculated  in  order  to  determine  the  CF.  Overall,  the  CF  depends 
on  the  correlation  of  the  inverse  of  the  initial  downrange  velocity 
with  the  initial  azimuth  velocity  and  with  the  crosswind.  Equa¬ 
tions  (53)  and  (54)  require  the  information  of  vx  ,  vyQ,  u,Q,  and  wx. 
In  order  to  better  estimate  their  correlations,  a  good  number  of 
experimental  data  are  necessary.  The  tediousness  and  complexity 
of  data  collection  for  all  three  velocity  components  and  time- 
dependent  wind  speeds  for  a  sizable  number  of  shots  could  hinder 


using  the  formulation.  In  some  cases,  estimation  from  statistical 
inference  may  be  necessary.  In  other  cases,  historical  data  for  a 
weapon  system  could  be  utilized  to  compensate  the  difficulty. 

No  complete  experimental  data  for  low-speed  projectiles  are 
available  for  demonstration  here.  On  one  hand,  the  experiment 
that  was  designed  to  capture  information  required  for  the  CF  cal¬ 
culations  was  very  limited.  Much  experimental  data  were  insuffi¬ 
cient  to  provide  a  thorough  study  about  this  matter.  On  the  other 
hand,  significant  portions  of  the  projectiles  in  this  category  are  the 
ammunition  in  which  the  dispersion  data  are  fairly  sensitive.  The 
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Table  6  List  of  possible  values  of  C''d  for  low-speed  projectiles 


Combination 

cd 

m 

P 

cd 

i 

2.7 

0.179 

1.204 

0.01178 

2 

2.1 

0.179 

1.204 

0.00916 

3 

3.3 

0.179 

1.204 

0.01440 

4 

2.7 

0.174 

1.204 

0.01212 

5 

2.7 

0.184 

1.204 

0.01146 

6 

2.7 

0.179 

1.146 

0.01121 

7 

2.7 

0.179 

1.342 

0.01313 

inclusion  of  such  data  could  escalate  the  permission  issue  for  pub¬ 
lic  release.  Thus,  readers  are  left  to  utilize  Eqs.  (53)  and  (54)  for 
their  practices.  Nevertheless,  it  should  be  straightforward  to  fol¬ 
low  the  description  as  outlined  for  high-speed  projectile  applica¬ 
tion.  Generally,  it  is  difficult  to  directly  derive  the  quantities  of 

VarKo^o]’  Var^ro^'  VarK%1]’  ™d  Vartw*%1]-  Instead’  one 

can  collect  some  data  from  v ..  ,  vv  ,  v-  ,  and  wr  to  estimate  their 

xo  V  zo  x 

individual  statistical  properties.  By  leveraging  the  assumption  of 
Gaussian  distribution,  these  quantities  can  be  obtained  through 
Monte  Carlo  simulations. 

Alternatively,  a  notional  low-speed  projectile  is  adopted  for  the 
sensitivity  study.  It  possesses  a  mass  of  0.179  kg,  a  reference  area 
of  0.0013  m2,  and  a  drag  coefficient  of  2.7.  The  average  initial 
downrange  velocity  of  the  projectile  is  74.3  m/s,  with  a  standard 
deviation  of  0.794  m/s.  The  values  of  drag  coefficient,  projectile 
mass,  and  air  density  are  perturbed  to  determine  a  possible  range 
of  C*d.  The  calculated  Cd  values  are  given  in  Table  6.  It  can  be 
seen  that  the  drag  coefficient  is  the  biggest  driver  that  results  in 
maximum  and  minimum  Cd  values  of  0.01440  and  0.00916,  re¬ 
spectively.  The  extreme  values,  along  with  the  average  Cd,  are 
used  for  dispersion  study. 

Based  on  Eqs.  (53)  and  (54),  the  calculation  of  CF  requires  the 
quantity  of  Varf/T1  uy()],  Var[iT2],  Var[iT' uX()],  and  Var[i/^w], 
which  are  computed  to  be  6.79  X  10-3,  1.49  X  10-11,  1.70  X  10-3, 
and  3.47  X  10-3,  respectively.  A  wind  speed  with  a  standard  de¬ 
viation  of  4.5  m/s  is  utilized.  Since  no  closed-form  solution  is 
available  for  the  variances,  these  values  are  obtained  through 
Monte  Carlo  simulation.  A  total  of  seven  cases  with  different  pa¬ 
rameter  values  are  created  and  given  in  Table  7.  The  first  three 
cases  are  to  assess  the  effect  of  Cd  on  CF.  The  next  two  cases  are 
to  investigate  the  influence  of  Varlffi^1  j^]  and  Var[jT)1i'Xo]  on  CF, 
and  the  last  two  cases  are  to  study  the  impact  of  Var[iF2]  and 

Vartrrjiv]. 

The  CFs  at  five  different  ranges,  i.e.,  20  m,  30  m,  40  m,  50  m, 
and  60  m,  are  calculated  against  the  baseline  dispersion  at  10  m. 
The  calculated  CFs  for  the  elevation  and  azimuth  directions  are 
provided  in  Tables  8  and  9,  respectively.  It  was  found  that  correc¬ 
tion  on  angular  dispersion  in  the  elevation  direction  may  not  be 
needed  in  many  cases.  The  changes  in  Cd  have  little  effect  on  the 
CF.  The  primary  driver  is  the  variation  in  the  initial  downrange,  as 
shown  in  case  6.  The  larger  the  variance,  the  smaller  the  CF.  As 
discussed  previously,  a  CF  of  1.0  was  obtained  when  v-o  is  con- 


Table  7  List  of  parameter  values  for  the  case  of  varying  vZq 


Case 

C* 

Var[v_1i'v  ] 

L  zo  W 

var[i/.2] 

\ar[v7lvx  ] 

L  Zq  *oj 

YarO^rvJ 

i 

0.01178 

6.79  X10“3 

1.49  X  10_n 

1.70  X  10“3 

3.47  X10“3 

2 

0.00916 

6.79  X10“3 

1.49  X  10-11 

1.70  X  10“3 

3.47  X10“3 

3 

0.01440 

6.79  X10“3 

1.49  X  10-11 

1.70  X  10“3 

3.47  X10“3 

4 

0.01178 

6.79  X10“2 

1.49  X  10-11 

1.70  X  10“2 

3.47  X10“3 

5 

0.01178 

6.79  X10“4 

1.49  X  10-11 

1.70  X  10“4 

3.47  X10“3 

6 

0.01178 

6.79  X10“4 

1.49  XI  0“10 

1.70  X  10“4 

3.47  X10“2 

7 

0.01178 

6.79  X10“2 

1.49  XI  0“12 

1.70  X  10“2 

3.47  X10“4 
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Table  8  Correction  factors  against  dispersion  at  10  m  in  eleva¬ 
tion  direction 


Case 

20  m 

30  m 

40  m 

50  m 

60  m 

i 

1.000 

1.000 

1.000 

1.000 

1.000 

2 

1.000 

1.000 

1.000 

1.000 

1.000 

3 

1.000 

1.000 

1.000 

1.000 

1.000 

4 

1.000 

1.000 

1.000 

1.000 

1.000 

5 

1.000 

1.000 

0.999 

0.998 

0.997 

6 

0.999 

0.996 

0.992 

0.985 

0.974 

7 

1.000 

1.000 

1.000 

1.000 

1.000 

Table  9  Correction  factors  against  dispersion  at  10 
muth  direction 

m  in  azi- 

Case 

20  m 

30  m 

40  m 

50  m 

60  m 

i 

0.988 

0.965 

0.933 

0.892 

0.842 

2 

0.993 

0.980 

0.961 

0.936 

0.905 

3 

0.981 

0.947 

0.898 

0.838 

0.770 

4 

0.999 

0.996 

0.993 

0.987 

0.980 

5 

0.899 

0.771 

0.647 

0.541 

0.455 

6 

0.639 

0.440 

0.325 

0.252 

0.203 

7 

1.000 

1.000 

0.999 

0.999 

0.998 

stant.  In  the  azimuth  direction. 

a  lower  value  of  Cd  improves  CF. 

At  the  range  of  60  m, 

a  high  drag  coefficient  may  result 

in  a  CF  as 

low  as 

0.77.  Overall, 

the  increase  in  the  variation  of  the  initial 

azimuth  velocity  has  a  limited  impact  on  CF.  The  major  driver  is 
the  variance  of  the  product  of  the  crosswind  velocity  and  the 
inverse  of  the  initial  downrange  velocity.  One  can  observe  from 
cases  5  and  6  that  when  Var[iT|  vw  ]  is  substantially  less  than 
VarfjrJj^],  the  CF  may  deteriorate  to  a  surprisingly  low  level  of 
0.203.  Since  the  increase  in  the  variation  in  initial  downrange 
velocity  magnifies  both  variances  of  the  products  simultaneously, 
a  marginal  impact  on  CF  is  expected.  As  a  result,  the  variation  in 
crosswind  velocity  appears  to  be  dominant,  which  is  in  line  with 
the  conclusion,  as  described  in  the  application  of  high-speed  pro¬ 
jectiles. 

6  Conclusion 

This  study  formulated  ballistic  dispersion  based  on  the  assump¬ 
tion  that  the  spread  of  TIPs  is  primarily  driven  by  the  variation  in 
the  initial  velocity  of  projectiles  as  they  exit  a  gun  barrel.  In  this 
study,  the  variation  in  the  initial  velocity  was  considered  to  be  a 
Gaussian  distribution  in  both  elevation  and  azimuth  directions. 
The  formulation  started  with  a  simplistic  case,  i.e.,  an  ideal  situ¬ 
ation  where  no  external  forces  are  imparted  to  the  projectiles  over 
the  entire  flight  period.  This  scenario  agrees  with  the  conventional 
calculation  of  ballistic  dispersion  in  angular  mil.  The  angular  dis¬ 
persion  was  obtained  from  dividing  the  standard  deviation  of  TIP 
by  the  corresponding  downrange  distance.  The  dimensionless  unit 
enables  us  to  make  dispersion  comparisons  at  multiple  downrange 
distances.  Simply  speaking,  the  comparison  was  based  on  a  linear 
relationship  between  the  absolute  dispersion  and  the  downrange 
distance,  and  the  relationship  was  illustrated. 

The  formulation  of  dispersion  was  extended  to  a  more  realistic 
situation  where  air  drag,  gravity  drop,  and  crosswind  were  taken 
into  account.  Projectiles,  of  which  the  initial  downrange  velocity 
may  be  considered  as  a  constant  from  one  round  to  another,  were 
first  investigated.  In  the  elevation  direction,  dispersion  was  found 
to  be  linearly  proportional  to  the  downrange  distance,  whose  re¬ 
sult  agrees  with  the  ideal  case.  In  other  words,  air  drag  and  gravity 
make  no  contribution  to  dispersion.  In  the  azimuth  direction,  the 
drift  distance  is  expected  to  be  zero.  However,  dispersion  is  al- 
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ways  magnified  because  of  crosswind  regardless  of  the  air  flow 
direction.  The  larger  the  variance  of  crosswind  velocity,  the 
greater  the  dispersion. 

Furthermore,  dispersion  was  formulated  for  projectiles  that  pos¬ 
sess  varying  initial  downrange  velocity  from  one  round  to  another, 
which  is  even  closer  to  reality.  The  variation  in  the  velocity  was 
modeled  by  a  Gaussian  random  variable,  the  same  representation 
as  the  initial  elevation  and  azimuth  velocities.  In  this  case,  the 
flight  time  to  a  given  downrange  distance  is  uncertain,  owing  to 
the  stochastic  property  of  the  projectile  launch  velocity.  As  a  re¬ 
sult,  the  variation,  in  conjunction  with  gravity  effect,  becomes  a 
significant  factor  in  determining  dispersion  in  the  elevation  direc¬ 
tion.  In  the  azimuth  direction,  three  random  variables  were  in¬ 
volved  in  the  dispersion  formulation.  Thus,  a  complex  expression 
was  derived,  which  requires  the  computation  of  the  correlations  of 
the  initial  downrange  velocity  with  the  initial  azimuth  velocity 
and  with  the  crosswind  velocity.  While  sizable  amounts  of  data  to 
determine  the  correlations  may  not  be  feasible,  the  statistics  of 
each  individual  factor  can  be  estimated  with  a  few  data  points, 
followed  by  comprehensive  Monte  Carlo  simulation  for  the  cor¬ 
relation  calculation. 

In  this  study,  variant  physical  conditions  were  utilized  to  for¬ 
mulate  general  ballistic  dispersion  from  a  stochastic  approach. 
The  theoretical  formulation,  which  links  initial  projectile  veloci¬ 
ties  to  target  impact  distribution,  may  be  useful  to  facilitate  the 
assessment  of  weapon  performance.  For  instance,  with  existing 
firing  data,  one  can  estimate  the  initial  elevation  and  azimuth  ve¬ 
locities  for  the  weapon  system.  On  the  other  hand,  when  the  pro¬ 
jectile  exit  conditions  at  the  muzzle  are  specified  with  a  weapon/ 
ammunition  system,  it  would  be  straightforward  to  determine 
dispersion  outcome  in  both  elevation  and  azimuth  directions. 

The  derived  formulations  confirm  a  highly  nonlinear  relation¬ 
ship  of  dispersion  with  the  downrange  distance  in  several  cases. 
While  the  traditional  conversion  to  dispersion  in  angular  mil  is 
quick  and  easy,  the  quantity  tends  to  be  overestimated,  particu¬ 
larly  at  longer  distances.  Thus,  CFs  are  proposed  for  dispersion 
comparison  at  multiple  downrange  distances.  Based  on  a  paramet¬ 
ric  study,  CFs  were  found  to  be  nontrivial  in  several  cases.  In  the 
elevation  direction,  the  variation  in  initial  downrange  velocity  was 
shown  to  be  the  primary  contributor  to  the  CF.  The  effects  of  drag 
coefficient,  projectile  mass,  and  air  density  were  very  limited.  In 
the  azimuth  direction,  the  variation  in  crosswind  velocity  was  de¬ 
termined  to  be  the  major  driver  for  the  CF.  The  drag  coefficient, 
projectile  mass,  and  air  density  were  found  to  slightly  influence 
the  CF.  Through  this  quantitative  study,  one  can  better  understand 
the  degree  of  impact  on  ballistic  dispersion  due  to  variant  projec¬ 
tile  conditions  at  the  muzzle  and  external  force  fields. 
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Nomenclature 

x  =  azimuth  direction 
y  —  elevation  direction 
z  =  downrange  direction 
Xq  =  gun  muzzle  position  at  firing  in  the  azimuth 
direction 


y  o 


Zo 


jp 

yj 

Zp 

dr 


Hr 


“y0 

v*o 

Vx 


8 

E[n] 

Varjrt] 

sd[ii] 
Cov[fIi  ,tl2] 


gun  muzzle  position  at  firing  in  the  elevation 
direction 

gun  muzzle  position  at  firing  in  the  downrange 
direction 

MPI  position  in  the  azimuth  direction 
MPI  position  in  the  elevation  direction 
MPI  position  in  the  downrange  direction 
separation  distance  of  an  impact  point  from 
MPI  in  the  azimuth  direction 
separation  distance  of  an  impact  point  from 
MPI  in  the  elevation  direction 
downrange  distance  from  the  firing  position 
(AoDO’Zo) 

separation  distance  between  the  firing  point 

and  the  aim  point  in  the  elevation  direction 

initial  velocity  of  a  projectile  in  the  azimuth 

direction  at  a  gun  muzzle 

initial  velocity  of  a  projectile  in  the  elevation 

direction  at  a  gun  muzzle 

initial  velocity  of  a  projectile  in  the  downrange 

direction  at  a  gun  muzzle 

flight  velocity  of  a  projectile  in  the  azimuth 

direction 

flight  velocity  of  a  projectile  in  the  elevation 
direction 

flight  velocity  of  a  projectile  in  the  downrange 
direction 

wind  velocity  in  the  azimuth  direction 
flight  time 

gravitational  acceleration 
expected  value  of  a  random  variable  D 
variance  of  a  random  variable  Cl 
standard  deviation  of  a  random  variable  Cl 
covariance  between  two  random  variables  11  [ 
and  Cl2 
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